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ABSTRACT 

Bright single and binary stars were observed at the 4.1-m telescope with a fast electron- 
multiplication camera in the regime of partial turbulence correction by the visible-light adaptive 
optics system. We compare the angular resolution achieved by simple averaging of AO-corrected 
images (long-exposure), selection and re-centering (shift-and-add or "lucky" imaging) and speckle 
interferometry. The effect of partial AO correction, vibrations, and image post-processing on the 
attained resolution is shown. Potential usefulness of these techniques is evaluated for reaching 
the diffraction limit in ground-based optical imaging. Measurements of 75 binary stars obtained 
during these tests are given and objects of special interest are discussed. We report tentative 
resolution of the astrometric companion to ( Aqr B. A concept of advanced high-resolution 
camera is outlined. 

Subject headings: Astronomical Instrumentation 
1. Introduction 



([Bates fc Cadv I I198QI : iBates fc Frightl Il982[ ) has 



Progress of observational astronomy depends 
mostly on advances in the wavelength coverage, 
sensitivity, and angular resolution. Although the 
record in resolution is set by long-baseline interfer- 
ometers, imaging and spectroscopy of faint sources 
is done today at the angular resolution of adaptive 
optics (AO) in the infra-red (IR) and at the res- 
olution of the 2.4-m Hubble Space Telescope in 
the visible. Reaching the diffraction limit in the 
visible from the ground for progressively fainter 
sources continues to be an important goal for as- 
tronomy. Here we report on new experiments in 
this direction. 

Diffraction-limited resolution at large ground- 
based telescopes has been f irst a ttained by speckle 
interferometry (Labeyriel Il970l). In the follow- 
ing years, modifications to this method to obtain 
true images, rather than their auto-correlations 
(ACFs), have been developed. One of the early 
approaches consisted in shifting the images to con- 
cidence of their brightest pixels and accumulat- 
ing the result. This shift-and-add (SAA) method 



been used to reconstruct the first image_of stel- 
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lar surface ( Lvnds. Warren, fc Harvev J 119761). It 
has been demonstrated that this appro ach is su- 
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erior to simple tip-tilt compensation (jChristou 



1991). Tip-tilt correction works well only at mod- 
erate D/ro ratios {D - telescope diameter, ro - 
Fried parameter); in this case, selection of the 
sharpest images with prominent central speckle is 
a good strategy. For D/ro > 10, the probablity of 
obtaining diffraction-limited image without turbu- 
lence correction becomes negligible, but the SAA 
method and speckle interferometry still work. 

The development of electron-multiplication 
CCD (EMCCD) detectors made it practical to 
record sequencies of short-exposure frames with- 
out readout noise penalty. On this occasion, the 
SAA met hod was revived under a new ndime, lucky 
imaging (iBaldwin. Warner, fc McCav I [20081 ). A 
combination of re-centering on bright pixel and se- 
lection of best frames was proposed. The selection 
works well only at small D/tq^ restricting this ap- 
proach to long wavelengths and to telescopes with 
D < 2.5 m. However, partial turbulence correc- 



tion with AO effectively increases ro, and the SAA 
method can give good results at larger apertures 
( Law et aI~ll2QQ9[ ). Dedicated high-resolution im- 
agers such as LuckyCam (iLaw et al. I [200^ or 
FastCam (|Qscoz et al. I l2QQ8[ ) were built to take 
advantage of these new opportunities. 

An AO system for partial image correction in 
the vsible, SOAR Adaptive Module (SAM ) is be- 
ing built at CTIO (IXokovinin et al. II2QQ41 l2QQ8al ) 
for the Southern Astrophysical Research Telescope 
(SOAR). This instrument will use a Rayleigh laser 
guide star to selectively correct low turbulence lay- 
ers. The diffraction limit will not be reached, 
but the effective increase of ro resulting from 
the partial correction will boost the potential of 
other methods such as speckle interferometry and 
SAA. First t ests of these me thods at SOAR are 
described bv iTokovmin fc Can tarutti ( 2008^) and 
Cecil fc Rashkeev I (|2QQ7[ ). respectively. 

In this paper we report some results from 
the first tests of SAM. At this stage, natu- 
ral guide stars (NGS) were used to close the 
AO loop. We begin with the description of 
our experiments in Section [2] and the analy- 
sis of the AO performance in Section [3l Then 
three methods of data processing - simple av- 
eraging, SAA, and speckle interferometry - are 
compared in Section HJ Measurements of close 
binary stars obtained as a by-product of these 
tests are reported in Section O they continue 
the series of speckle observati on s publ ished by 
Tokovinin. Mason, fc Hartkopf I (|2QlQal ) (here- 
after TMHIO). Some pairs deserve special com- 
ments; two faint tertiary companions in wider 
binaries were resolved for the first time. We con- 
clude in Section [6] with a perspective of diffraction- 
limited imaging in the visible. 

2. Description of the experiments 

SOAR is located at the Cerro Pachon observa- 
tory in Chile, at 2800 m above the sea level. It has 
alt-azimuth mount with two Nasmyth and several 
bent Cassegrain focal stations fed by a flat mirror 
M3. This mirror can be actuated in tip and tilt. 

The SOAR Adaptive Module (SAM) is an AO 
system for correcting wave-front distortions and 
feeding improved imag e s to s cience instruments 



(Tokovinin et al. 200 



age s to s cK 
3 l2008al) . 



SAM has inter- 
nal re-imaging optics (two off-axis paraboloids). 



turbulence corrector (60-element curvature de- 
formable mirror, DM), and a 10x10 Shack- 
Hartmann wave- front sensor (WFS). The SAM 
instrument is designed to work with a laser guide 
star, but initially it used natural guide stars 
(NGS). The light of an on-axis star was mostly 
reflected to the WFS by a neutral beamsplit- 
ter placed after the DM, with 14% transmitted 
to the science detector. The latter was a High- 
Resolution Camera (HRCAM) which contained 
an EM CCD with a pixel scale of 15.23 mas and 
filters (|Tokovinin fc Cantaruttn [2008). Owing to 
the 2.1^ light attenuation by the beamsplitter, we 
observed only relatively bright stars. 

SAM was installed in August 2009 at the Nas- 
myth focus, receiving the light after four reflec- 
tions in the telescope. It rotates as a whole for 
compensating field rotation as the telescope tracks 
stars. The AO loop cycle was 4.3 ms, with a 2- 
cycle delay between WFS measurement and cor- 
rection. To compensate for this relatively slow 
loop, we used the Smith predictor (SP) controller 
(Madec I| l998l ). Alternatively, a standard leaky in- 
tegrator (INT) with a gain of 0.25 was tried. The 
degree of AO correction is equivalent to compan- 
sation of 45 Zernike modes. 

Of the three engineering runs assigned to these 
tests in August - October 2009, the first was 
completely lost to clouds, the second produced 
1.5 nights of data under fast and poor seeing of 
1'' — 1.6'^ the third and last run enjoyed few hours 
of good conditions (seeing 0^5 to 1^'). 

We pointed telescope to single or multiple 
bright stars, closed the AO loop and recorded, 
typically, 400 images with 20-ms exposure time 
and EM gain of 44. Only the central region of 
3^' X 3^' was recorded in most cases. In parallel, 
we frequently saved data from the AO loop (se- 
quences of centroids and DM voltages) for off-line 
analysis. 

3. Analysis of the AO loop data 

We calculate the first 45 Zernike coefficients aj 
of the input (atmospheric) and corrected (resid- 
ual) wave- fronts (in radians at A = 0.5/im) from 
the AO loop data. The DM voltages are converted 
into corresponding slopes by matrix multiplica- 
tion, using the recorded AO interaction matrix. 
The slopes, in turn, are converted into Zernike 
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time To is derived from the temporal variation 
of the defocus coefficient a4, as des cr ibed b y 



Tokovinin. Keherer. Coude du Forestol (l2QQ8bh . 



Its values at A = 0.5/im ranged from 1 ms to 
3.5 ms. Both tq and tq derived from the loop 
data are in reasonable agreement with simulta- 
neous measurements from the MASS-DIMM site 



monitor (jTokovinin et al. ll2QlQbl ). 
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Fig. 1. — Temporal power spectra of Zernike co- 
efficients a4 (defocus, top) and ag (coma, bottom) 
recorded on October 2, 2009 under good seeing 
(ro = 0.13 m) in closed loop with the SP controller. 
Dash-dotted straight lines indicate the power-law 
slopes k = 2.3 and k = 3.0, respectively. The spec- 
trum of a4 has lines at 47 Hz, 50 Hz, and 94 Hz. 



coefficients by means of the calculated theoreti- 
cal slopes of Zernike modes (gradient matrix); we 
simply multiply the vector of slopes by the pseudo- 
inverse gradient matrix. Of course, the Zernike 
coefficients aj derived in this way contain errors 
due to noise on the measured slopes and cross-talk 
with higher-order modes (aliasing). 

The Fried radius ro and seeing sq = 0.98A/ro 
are estimated from the variance of Zernike co- 



efficients (j| = (Aa^) which should, in theory. 



equal Cjj {D/r^f^^ (INoII I Il976[) . The ratio of 



measured variances to the Noll's coefficients Cjj 
is averaged using modes from 7 to 19 to get an 
estimate of D/tq. The atmospheric coherence 



Fig. 2. — The 50-Hz vibration in tip and tilt. The 
signal from 48 Hz to 52 Hz is isolated by the tem- 
poral bandpass filter applied to the a2 and as co- 
efficients, converted to angular units, and plotted 
as XY trajectory of the optical axis. 

Temporal power spectra of the Zernike coeffi- 
cients corresponding to the input (un-corrected) 
wave-fronts are calculated from the DM signal, 
with division of the raw power spectrum by the 
theoretical frequency response of the closed AO 
loop. The high-frequency part of those spectra 
shows a typical power-law decline proportional to 
with exponents k between 2 and 3 (Fig. [T]). 
A high-frequency asymptotic decline with k = 
11/3 is characteristic of atmospheric tilts, while 
a steeper dependence with k = 1 7/3 is expected 



for the defocus and higher modes (jRoddier et al 
[1993). This mismatch between theory and real 
spectra has also been observed with physically 
simulated turbulence in the laboratory and repro- 
duced in the numerical simulations. Most proba- 
bly, it is caused by t he atmospheric errors of WFS 
slope measurements (Thomas et al. |[2'006[ ) and by 
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Fig. 3. — Comparison of various processing algorithms for the case of the binary star BU 172 (0'/408, 
V = 6.4; 6.6). The left column shows average images, the middle column - SAA images with 20% selection, 
the right column - ACF obtained in the speckle processing of the same data cube. The rows show three 
data cubes obtained one after another in a rapid sequence. Top row - / filter, INT controller; middle row 
- Ha filter, INT controller; bottom row - Ha filter, SP controller. The insert in the middle row shows the 
50-IIz vibration tajectory from the simultaneous loop data, on a different scale (cf. Fig. [2]) which matches 
the elongation of the ACF peaks in angle and amplitude. The PSF in the upper row is elongated by the 
atmospheric dispersion. 



the aliasing of higher-or der modes, as shown by 



Veran fc Poyneer 



(|2010h . 

The power spectra of the residual coefficients 
derived from the WFS signal show that distor- 
tions with low temporal frequencies (up to 20 Hz) 
are compensated, while high frequencies are am- 
plified. This amplification is caused by the AO 
servo system and depends on the type of temporal 
controller used; the SP controller gives more high- 
fequency errors than INT. The spectra of residual 
coefficients match very well their expected shape 
- a product of the input disturbance and the AO 
loop error transfer function (Fig. [T]). 

Quasi-periodic signals caused by vibrations are 
seen as narrow "lines" superposed on the other- 
wise smooth spectra of input or residual distor- 
tions. Such lines are detected in modes up to 6 
(astigmatism), but are not present in the higher- 



order modes. The focus coefficient almost al- 
ways has a line at 47 Hz. This frequency is too 
high to be caused by a real focus variation in the 
telescope. It is suspiciousy close to the 1/5 of the 
loop frequency and is likely of instrumental origin 
(we do not see such signal in the laboratory, how- 
ever). On the other hand, the 50-Hz component 
in the tip and tilt coefficients is real. This vibra- 
tion is caused by the pick-up noise in the position 
sensors of the SOAR fast tip-tilt mirror, M3. This 
signal is detected by the standard SOAR guider, 
as well as directly in the voltages of the M3 drives 
(M. Warner, private communication). By extract- 
ing the 50-Hz component from the SAM loop sig- 
nal, we show that motions on both axes are co- 
phased (Fig. [2]). The amplitude of these vibrations 
and relative phase shift between the axes are not 
constant. On October 2, 2009, the worst-case rms 
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amplitude was 19 mas, the smallest one - 3 mas. 

As the typical exposure time in speckle inter- 
ferometry is 20 ms - one full period of the 50- 
Hz vibration - the speckle structure is sometimes 
blurred and the resolution does not reach the 
diffraction limit (30 mas in the visible). If the vi- 
bration ellipse is narrow, as in Fig. [21 the resulting 
blur can mimic a close binary with separation of 
20-50 mas. Such effect was indeed seen in the pre- 
vious speckle data obtained at SOAR (TMHIO). 
It biases the measurements of very close binaries 
in un-controlled way. Modification of the M3 drive 
electronics is in progress. Once the 50-Hz compo- 
nent is eliminated, there will be no other sources of 
fast blur, making the SOAR telescope most suit- 
able for diffraction-limited observations. 

4. Imaging 

Series of short-exposure images were processed 
in three different ways (Fig. [3]). In all cases we 
subtracted the average background to remove the 
fixed offset and the dark current which was non- 
negligible in a small number of "hot" pixes. All 
images are Nyquist-sampled. 

Speckle interferometry in its standard form 
is implemented as described in TMHIO. We com- 
pute the power spectrum from all images in each 
series, remove the photon- noise bias, filter out 
the low frequencies and derive by inverse Fourier 
transform the average auto-correlation function 
(ACF). The phase information is lost, but the 
diffraction-limited resolution can be reached even 
without AO correction. 

We found that fast wave-front distortions am- 
plified by the SP controller in closed loop (cf. 
Fig. [1]) blur the high-frequency structure in the 
20-ms images and degrade the resolution in the 
speckle mode. At the same time this controller 
gives similar or smaller overall wave-front residu- 
als as the INT controller with a gain of 0.25; the 
resolution in the long-exposure images is also a 
little better with the SP controller. 

All data on binary stars were processed uni- 
formly by speckle interferometry to derive rela- 
tive component positions and magnitude differ- 
ences (see Sect. [5]). For the analysis of AO imaging 
performance, we de-convolve the binary signature 
(as determined by the speckle processing) and de- 
rive the point spread function (PSF). 



Long-exposure images were obtained as sim- 
ple averages of the data cubes. Typically, the to- 
tal exposure was about 12 s (30 Hz frame rate, 400 
frames). As an option, we could re-center indi- 
vidual images before averaging; the images were 
shifted by integer number of pixels as determined 
by their centroids. Re-centering brings almost no 
resolution gain when the AO loop is closed. There- 
fore we restrict further analysis to long-exposure 
images. 

SAA imaging is implemented by selecting a 
certain fraction of the sharpest images in a given 
series and co-adding them with re-centering on the 
brightest pixel. In this case, the most prominent, 
diffraction-limited speckle is enhanced in the av- 
erage image, while fainter speckles are averaged 
out and contribute to the diffuse halo. As a re- 
sult, a PSF with a diffraction-limited core and a 
halo is obtained, reminiscent of partially-corrected 
AO images. Combination of partial AO correction 
with re-centering on the brightest pixel gives bet- 
ter results than each of these methods alone. 

For the image selection, we experimented with 
different criteria such as second moment or sharp- 
ness and found that selection on the maximum 
image i ntens i ty wo rks best, in agreement with 
[Smith et al. ( 20091 ). As to the fraction of retained 
images, it is not critical. The resolution gain is 
achieved mostly by re-centering, not by selection. 
We co-add the best 20% frames to produce the 
SAA images. For faint stars, the choice of the 
brightest pixel is affected by the noise, leading 
to a sharp 1-pixel spike at the center; this spike 
was replaced by the average of its neighbours for 
PSF evaluation. In the case of a binary star with 
components of comparable brightness, the noise 
and speckle fluctuations cause the brightest pixel 
to belong sometimes to the secondary component. 
As a result, a weaker false spike located symmet- 
rically with respect to the primary component ap- 
pears in the SAA image. 

Figure |4] gives some representative images of 
multiple systems observed on October 2, 2009 
in closed loop with the SP controller. As ex- 
pected, the best resolution in the long-exposure 
compensated images is reached in the / band 
(wavelength 0.8 /im), so most of the data were 
obtained in this filter. The resolution gain over 
un-compensated, seeing-limited imaging is impres- 
sive, even at shorter wavelengths. For example. 
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FIN 329 0.12" I 385 0.27"+0.40" HR 695 0.50" 

Fig. 4. — Images of some multiple systems in the 
/ band. Left: a bright close pair FIN 329, long- 
exposure. Center: tight visual triple system I 385 
with comparable separation between components, 
long exposure. Right: astrometric binary Hi For 
with a faint (Am = 3.4) companion first resolved 
at SOAR in 2007 (SAA image). 



three images taken one after another in the /, V, 
and B filters show the long-exposure FWHMs of 
0.16'', 0.30'' and 0.44" respectively under a 0.80" 
seeing. 



better resolution than simple averaging. However, 
as the seeing improves and the AO correction be- 
comes good, the additional gain brought by this 
technique becomes less evident (Fig. [5j). Average 
FWHM for the 22 points in Fig. [5] is 0'.'25 and 
0'.'18 for the long-exposure and SAA, respectively. 
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Fig. 5. — Relation between atmospheric seeing 
(horizontal axis) and the FWHM of the compen- 
sated images in the case of long exposure (squares) 
and SAA (crosses connected to squares) process- 
ing. Data of October 2, 2009, /-band. 

The SAA (co-addition of re-centered images) 
produces a FWHM with diffraction-limited core 
even without AO correction. We obtained such 
images from the data of 2008-2009 speckle runs 
at SOAR recorded under a O'.'S seeing. Par- 
tial AO correction helps to concentrale the light 
in few bright speckles and thus enhances the 
range of condition s where this method works well 




( Law et al. I l2009h . We found that under partial 



AO correction the SAA post-processing yields a 
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Fig. 6. — Radial profiles of the PSFs of long- 
exposure (full lines) and SAA (dashed lines) im- 
ages under good (top, ro = 0.18 m) and mediocre 
(bottom, ro = 0.095 m) seeing in the / band. 
Gaussian PSFs of seeing-limited width are plot- 
ted for reference. 

Radial profiles of the PSFs under good and par- 
tial AO correction are plotted in Fig. [6l For com- 
parison, Gaussian profiles with the seeing-limited 
FWHM 0.98A/ro are given, using ro values esti- 
mated form the AO loop. The shape of the AO- 
corrected PSFs is different from a Gaussian and 
resembles a negative exponent. Such PSF is typi- 
cal for an AO system like SAM, according to our 
simulations. 

A detailed comparison of various processing 
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techniques is presented in Fig. [3l The binary star 
BU 172 was observed at zenith distance 29.5° un- 
der average seeing (ro = 0.12m to 0.14m from the 
loop data). The first row demonstrates good com- 
pensation in the / band (long-exposure FHWM 
O'.'IS) achieved with the INT loop controller. The 
SAA and speckle processing increase the resolu- 
tion to the point where the atmospheric dispersion 
(47 mas in the direction parallel to the binary) be- 
comes visible. The dispersion is negligibly small 
(5 mas) in the narrow-band Ha filter, but the AO 
compensation at this longer wavelength is not so 
good (FWHM 0V33). The photon flux is low, so 
the SAA method does not work very well, while 
the speckle processing gives a stronger signal. The 
images are elongated perpendicularly to the bi- 
nary by the 50-Hz vibrations (see the insert) of 
peak-to-peak amplitude 45 mas, in this case. Fi- 
nally, when the SP controller is used, the images 
become more fuzzy because of the rapid residual 
wavefront jitter. However, the average image is 
even a little sharper (FWHM 0'.'31) than with the 
INT controller. 

Although the binary BU 172 is bright, the flux 
is decreased by 2.1^ by the beamsplitter in SAM 
and, further, by 3.4^ by the narrow-band Ha filter 
(in comparison with V or R bands). Therefore, 
the photon flux is equivalent to a binary star with 
V = 11.9,12.1 observed in a wide band without 
SAM beamsplitter. Further gain in sensitivity by 
at least ~ 1^ is achievable with a better detector; 
the quantum efficiency is presently only about 0.5, 
while the clock injection charge (about 40 photons 
per frame) is a serious limitation for faint stars. 
Considering that this observation is not yet close 
to the limit, there are reasons to believe that stars 
down to V ~ 15^ can be observed at diffraction- 
limited resolution in speckle or lucky mode with 
SAM. 

5. Results on multiple stars 
5.1. Binary-star measurements 

Relative positions and magnitude difference in 
close binary stars were measured by speckle in- 
terferometry, following the same procedure as in 
TMHIO. We took advantage of that work to cali- 
brate new results from SAM on common wide sys- 
tems and found that the pixel scale remained the 
same within errors, 15.23 mas. An offset in posi- 



tion angle of 6.0° has been determined and cor- 
rected. For the most part, we measured bright 
and well-known systems, although some recently 
discovered close pairs were also observed. 

Table [1] lists 122 measurements of 75 resolved 
binary stars. Its co lumns contain (1) the WDS 



(|Mason et al. Il200ll ) designation, (2) the "discov- 
erer designation" as adopted in the WDS, (3) an 
alternative name, mostly from the Hipparcos cat- 
alog, (4) Besselian epoch of observation, (5) filter, 
(6) number of individual data cubes, (7,8) position 
angle in degrees and internal measurement error 
in tangential direction pag in mas, (9,10) separa- 
tion p and its internal error ap in mas, and (11) 
magnitude difference Am. An asterisk follows if 
Am is determined from the resolved photometry 
(see TMHIO for details). A colon indicates that 
the data are noisy and Am is likely over-estimated. 
Note that in the cases of multiple stars, the po- 
sitions and photometry refer to the pairings be- 
tween individual stars, not with photo-centers of 
sub-systems. 

For stars with known orbital elements, columns 
(12-14) of Table □ list the residuals to the 
ephemeris position and th e reference to the orbi t 
from the Orbit Catalog ([Hartkopf et al. Il200l[ ). 

Some measurements from t his work were used to 

calcul ate improved orbits bv lMason. Hartkopf. fc Tokovinin 
(|201 0V these orbits are referenced as Msn2010 in 
Table E 

In most cases we see binary companions in the 
SAA images and can thus resolve the 180° ambi- 
guity inherent to the speckle method. For a few 
pairs where this was not possible, the position an- 
gles are preceded by the letter A (ambiguous). 

5.2. Comments on individual systems 

00522-2237 = STN 3 has a new faint tertiary 
component D at 0^57 from B with Al = 4"^. As 
the separation in the wide pair AB is only 2'', the 
triple system appears to be close to the dynamical 
stability limit. We cannot exclude that the new 
companion is just an unrelated background source 
projecting on the AB. The companion C listed in 
the WDS is in fact such a background source, as 
evidenced by its relative motion since 1897. How- 
ever, the actual separation of AC is > 10'', so the 
C companion was well outside the field. 

02057-2423 = HIP 9744 is a close visual 
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Fig. 7. — Preliminary visual orbit of HIP 9774 
Aa,Ab = WSI 71. 



triple star consisting of the 0^'8 pair AB = 
I 454 and the inner spectros copic binary Aa,Ab 
( Tokovinin fc Smekhov |[20Q2 ). It forms a quadru- 
ple or quintuple system together with HIP 9769 
which has common radial velocity, parallax, and 
proper motion with HIP 9774, at 56^' separation. 

We determined a preliminary spectroscopic or- 
bit of the Aa,Ab sub-system with a period of 
2.6 y and eccentricity 0.76. The estimated semi- 
major axis of this pair is 43 mas, indicating that 
it can be resolvable at 4-m telescopes. First reso- 
lution of Aa,Ab (named on this occasion WSI 71) 
was achieved in July 2008, shortly after passage 
through the periastron (TMHIO). The pair was 
resolved again in October 2008 at SOAR and ob- 
served also on October 2, 2009 with SAM. By that 
time the components opened up to 42 mas. 

An attempt to combine the three position mea- 
surements with the known spectroscopic orbit 
has shown that the separations and magnitude 
differences in the published speckle observations 
are likely over-estimated. They were determined 
jointly with the position of the wide pair Aa,B 
by fitting a triple-star model. We re-processed all 
speckle data by imposing a constraint Am = 0.30, 
as determined from the relative intensity of com- 
ponent's lines in the resolved spectra. The pres- 



ence of the wide companion B was ignored by fit- 
ting only the close binary. This reprocessing gives 
closer separations of Aa,Ab, as expected (Table [2]). 
No formal errors are given because measurements 
of such close binaries can be seriously biased by 
vibrations (see above). 

FigureElshows a preliminary interferometric or- 
bit obtained by fitting jointly radial velocities and 
resolved measurements. Clearly, more data are 
needed for a definitive solution. However, plau- 
sible values of the visual elements are obtained: 
semi-major axis a = 58 mas, position angle of 
nodes Q = 125° and inclination i = 77°. 

Speckle interferometry is the only method for 
resolving such close binaries. Adaptive optics at 
8-m telescopes has somewhat poorer resolution 
(diffraction limit X/D = 32 mas in the J band), 
whereas long-baseline interferometers lack the sen- 
sitivity and productivity to follow such binaries. 

22116-3428 = HIP 109561 is a visual triple 
consisting of the 0V4 pair AB = BU 769 and the 
inner pair Aa,Ab = CHR 230. We marginally re- 
solved the inner sub-system at almost the same 
position angle as measured in TMHIO (its esti- 
mated period is 100 y). Note that component A is 
a KlIII giant and B is a blue star, possibly of A5V 
spectral type. Therefore the magnitude difference 
of AB is larger at longer wavelengths. 

22288-0001 = C Aqr B = HR 8558 is a sec- 
ondary in the wide visual pair AB = STF 2909 
with components of similar spectral types and 
magnitudes (F3V and F6IV, V = 4.42; 4.51) 
observed since 1821. Both stars are above the 
Main Sequence and have fast axial rotation. A 
"wavy" motion of AB sugg e sted t hat there is an 
inner sub-system. Strand (|l942[ ) published the 
first astrometric orbit of this sub-system, attribut- 
ing it to the component B. Several authors con- 
firmed subsequently the 25-y astrometric wave and 
refined both outer and inner orbits; the latest 



Table 2: Measurements of HIP 9774 Aa,Ab 



Epoch 





P 


[0-C], 


[0-C], 




(deg) 


n 


(deg) 


n 


2008.5406 


132.4 


0.027 


-9.9 


0.002 


2008.7700 


188.7 


0.019 


15.8 


0.000 


2009.7558 


286.3 


0.042 


5.1 


-0.007 
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Fig. 8. — Fragments of the ACFs of ( Aqr B in the 
/ (left) and Ha (right) bands. The arrows mark 
the tentatively resolved astrometric companion Bb 
at 0^4. 



0.6 



(D 
(f) 
O 



0.4 



-0.0 



-0.2 



-0.4 



-0.6 




-0.6 -0.4 



-0.2 -0.0 0.2 

West, arcsec 



0.4 0.6 



Fig. 9. — Observed position of ( Aqr Ba,Bb 
(square ) in co mparison with astrometric orbits of 
Heintz I ([19841 . fuh line) and lScardia et al. I (|2QQ9l 
dashed Hne). For both astrometric orbits, the po- 
sition angle is inverted and the semi-major axis 
a =0^34 is adopted, as appropriate for direct res- 
olution. The companion is connected to its pre- 
dicted positions on both orbits. 



are pa pers bv lHeint"z" ( 1984 ) and Scardia et al. 
(120091). The first direct resolution of the astro- 



metric co mpanion in the K and L bands was re- 
ported bv lMcCarthv et al. I (|l982[ ) who fitted one- 
dimensional scans in the North-South direction to 
a triple-star model. No further resolved observa- 
tions of Ba,Bb were published to date. The com- 
ponent A was also resolved with speckle into a 
0V064 binary by lEbersberge r fc Weigelt I (Il979h . 
but never confirmed ever since, despite several 



speckle meas urements of A B at 4-m telescopes. 
Interestingly, iHeintz I ('1984') dismisses both resolu- 
tions as bogus and attributes the 25-y sub-system 
to A, based on the astrometry of A and B. 

Despite the controversy, the existence of a 25-y 
sub-system makes little doubt. Let us assume the 
mass of the Ba star Mpa = 1.5 in agreement 
with its spectral type and the dynamical mass sum 
derived from the AB orbit and adopt the mass sum 
of 1.9 for Bab. Then its period 25.8 y and the 
third Kepler's law lead to a semi-major axis of 
10.8 AU, or a =0^34 with the Hipparcos parallax 
of 31.5 mas. The comparison with the astromet- 
ric semi-major axis a = 75.5 mas (|Heintz 1984) 
leads to a mass ratio q = Meb/^Ba = 0.28 and 
^Bb = 0.4 M0, in broad agreement with other 
published estimates. The unseen component could 
be an MOV dwarf. In this case the magnitude dif- 
ference with Ba in the / band would be around 
rpj^^ estimated radial- velocity variation of Ba 
is Ki ~ 1.3km/s; it has never been detected, pos- 
sibly because of the fast axial rotation. 

We see a faint (Am ~ 5^) companion to B in 
four data cubes - two in the band / and two in 
Ha (Fig. [8|). This resolution should be considered 
as tentative until confirmed by further measure- 
ments. It would be relatively easy to detect such 
a companion with AO in the infrared. The com- 
ponent A has also been observed with SAM and 
found unresolved. 

The observed position of Ba,Bb, if it is real, 
does not match well the two recent astrometric 
orbits, which also differ from each other substan- 
tially (Fig. [9]). In both cases we adopted the semi- 
major axis of 0^34. If we correct the periastron 
epoch in the orbit of Heintz (|l984l ) from the orig- 



inal 2007.0 to 2004.0, it would agree with our ob- 
servation. 

6. Conclusions and outlook 

To summarize, the first tests of SAM in its NGS 
mode have shown the following: 

• The AO system works as expected when 
closed on natural guide stars. The FWHM 
resolution of AO-corrected long-exposure 
images can be better than 0'.'2 in the /-band 
and is significantly improved with respect to 
the seeing across the visible spectrum. 
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• Fine speckle structure in the images with 20- 
ms exposure time can be blurred by the 50- 
Hz vibration originating in the SOAR tip- 
tilt mirror (soon to be improved) and by the 
amplification of fast atmospheric distortions 
in the AO loop using the SP controller. A 
softer controller is preferable for diffraction- 
limited imaging. 

• The SAA ("lucky") processing of the data 
cubes improves the resolution for bright 
sources under partial compensation. It gives 
little additional gain under good AO correc- 
tion and is less powerful than the speckle 
processing in the photon-starved regime. 

• Both speckle and SAA work well for stars 
as faint as V ~ 12^ if all available photons 
are used by the current HRCam. A further 
gain of 1-2 magnitudes could be achievable 
by using state-of-the-art EM CCD with a 
higher quantum efficiency and lower charge 
injection noise. 

These first results suffer from small statistics 
and incomplete sampling of atmospheric condi- 
tions; more experiments are needed to substan- 
tiate and extend these findings. 

Methods like speckle and SAA imaging, en- 
hanced by partial AO correction, open a unique 
high-resolution niche. A non-exhaustive list of po- 
tential science programs includes: 

• Followup of close nearby binaries and de- 
termination of their orbits, leading to mass 
measurement (including pre-main-sequence 
stars), statistics of orbital elements, dynam- 
ics of multiple systems, etc. 

• Surveys of various stellar polulations for bi- 
narity delivering observational constraints 
on star formation. 

• Optical imaging of solar-system objects 
(planetary satellites, as teroids), see e.g. 
(jCecil fc Rashkeev II2007I ). 



Further interesting results are expected from 
the combination of the HRCam with SAM. In 
a more distant perspective, a dedicated high- 
resolution instrument will be worth considering. 
As the AO technology becomes mature and ac- 
cessible, we envision building a small-field camera 
with integrated compact AO system. The light 
from the object (or its bright stellar component) 
will be divided in wavelength between the science 
channel and the WFS, e.g. by sensing the beam re- 
flected from the science filter. An advanced WFS 
with a EM CCD would be able to operate on much 
fainter stars than is possible with SAM now. In 
fact, even the faintest stars will be useful. When 
the number of photons is not sufficient for fast 
turbulence correction, the AO system will work in 
the slow mode to compensate residual telescope 
aberrations (mostly defocus and tilts), delivering 
to the speckle camera truly seeing-limited images. 

We thank the SOAR telescope team for their 
essential help with the installation of the SAM 
instrument and observations. This research has 
made use of the CDS and SIMBAD services and of 
the Washington Double Star Catalog maintained 
at the U.S. Naval Observatory. 
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Table 1 

Measurements of binary stars at SOAR 



WDS 
(2000) 


Discoverer 
Designation 


Other 
name 


Epoch 
+2000 


Filt 


N 


e 

(deg) 


(mas) 


P 

n 


(jp 
(mas) 


Am [O-C]^ 
(mag) (deg) 


[o-c]p 

(") 


Reference 
code* 


00143-2732 


HDS 33 


HIP 


1144 


9 


6709 


I 


2 


310.8 


0.3 


0.1674 


0.9 


0.9 








00284-2020 


B 1909 


HIP 


2237 


9 


6709 




3 


307.8 


0.5 


0.1888 


0.5 


0.2 


-3.8 


0.003 


Sodl999 


00315—6257 


I 260 CD 


HIP 


2487 


9 


6711 




2 


242.4 


0.4 


0.2529 


0.3 


1.5 


5.6 


—0.017 


Msn2001c 


00345-0433 


D 2 AB 


HIP 


2713 


9 


6708 


I 


2 


A207.9 


0.4 


0.0610 


0.1 


0.8 


109.9 


-0.223 


Sey2002 










9 


6708 


R 


2 


A207.8 


1.0 


0.0604 


0.9 


0.1 


109.8 


-0.223 


Sey2002 










9 


7556 




6 


A204.7 


5.3 


0.0574 


1.2 


0.9 


106.6 


-0.228 


Sey2002 


00352-0336 


HO 212 AB 


HIP 


2762 


9 


6708 


I 


2 


229.6 


0.1 


0.2345 


0.4 


1.2 


-0.7 


0.000 


Msn2005 










9 


6708 




2 


229.5 


0.4 


0.2344 


1.0 


1.4 


-0.7 


0.000 


Msn2005 










9 


7556 




4 


231.7 


0.5 


0.2387 


0.2 


1.1 


-0.4 


0.000 


Msn2005 


00373-2446 


BU 395 


HIP 


2941 


9 


6709 




4 


95.0 


0.1 


0.4627 


0.1 


0.4 


2.8 


0.017 


Pbx2000b 










9 


7557 


I 


2 


95.6 


0.1 


0.4727 


0.3 


0.2 


* 2.9 


0.018 


Pbx2000b 


00427-6537 


I 440 


HIP 


3351 


9 


6711 




2 


269.4 


0.2 


0.4023 


0.3 


0.7 


1.1 


-0.084 


Lin2004a 


00522-2237 


STN 3 AB 


HIP 


4072 


9 


6709 




3 


243.1 


1.4 


1.9853 


2.3 


0.7 








00522-2237 


STN 3 BD 


HIP 


4072 


9 


6709 




3 


336.9 


5.9 


0.5743 


3.3 


4.0 








01061-4643 


SLR 1 AB 


HIP 


5165 


9 


6711 




2 


113.0 


0.1 


0.4050 


0.1 


0.5 


-29.2 


-0.042 


Ary2001a 


01078-4129 


RST 3352 


HIP 


5300 


9 


6711 




2 


156.1 


0.2 


0.1412 


0.3 


1.0 


1.9 


-0.003 


Sodl999 










9 


6711 




2 


156.1 


0.4 


0.1408 


0.3 


1.3 


2.0 


-0.004 


Sodl999 


01084-5515 


RST 1205 AB 


HIP 


5348 


9 


6711 




2 


110.0 


1.5 


0.5486 


0.4 


2.7 


-2.3 


0.033 


Lin2004a 


01094-5636 


HU 1342 


HIP 


5428 


9 


6711 




3 


328.8 


0.2 


0.3746 


0.7 


0.6 


-1.8 


-0.018 


Heil984a 


01158-6853 


I 27 CD 


HIP 


5842 


9 


6711 




2 


312.5 


0.3 


1.0540 


0.2 


0.5 


* 2.0 


-0.034 


Sodl999 


01198-0031 


FIN 337 BC 


HIP 


6226 


9 


6710 




2 


32.2 


1.4 


0.1337 


1.7 


0.8 


-3.8 


0.003 


Msn2010 


01198-0031 


STF 113 AB 


HIP 


6226 


9 


6710 




2 


19.1 


3.4 


1.6783 


0.3 


2.4 








01220—6943 


I 263 


HIP 


6377 


9 


6711 




2 


279.1 


0.2 


0.4271 


0.2 


0.8 


5.3 


-0.099 


Msnl999a 


01243-0655 


BU 1163 


HIP 


6564 


9 


6710 




3 


218.0 


0.2 


0.3227 


0.1 


0.4 


—0.4 


-0.002 


Sodl999 


01334-4354 


HDS 205 


HIP 


7254 


9 


6711 




2 


195.8 


0.2 


0.1322 


0.7 


0.9 








01350-2955 


BU 1000 AC 


HIP 


7372 


9 


6709 




2 


351.1 


2.1 


1.8128 


0.9 


3.2 


-17.8 


0.237 


Nwbl969a 


01350-2955 


DAW 31 AB 


HIP 


7372 


9 


6709 




2 


40.7 


2.6 


0.1366 


3.3 


0.5 


-7.0 


-0.005 


Msnl999c 


01361-2954 


HJ 3447 


HIP 


7463 


9 


6709 


I 


2 


182.5 


0.3 


0.7729 


0.1 


1.4 


-2.9 


-0.038 


Cve2006e 


01376-0924 


KUI 7 


HIP 


7580 


9 


6710 




1 


112.5 


0.2 


0.0771 


0.2 


1.2 


3.7 


0.002 


Tokl993 










9 


6710 




2 


111.9 


0.4 


0.0786 


0.7 


1.0 


3.1 


0.004 


Tokl993 


01417-1119 


STF 147 


HIP 


7916 


9 


6709 




2 


120.2 


0.1 


0.1072 


0.2 


1.0 
















9 


6710 




2 


120.4 


0.1 


0.1065 


0.5 


1.1 








02022-2402 


HDS 272 AB 


HIP 


9497 


9 


6710 




3 


340.6 


0.8 


0.5861 


1.6 


3.4 
















9 


6710 


R 


2 


340.6 


0.8 


0.5801 


0.3 


4.4 
















9 


7557 




3 


340.2 


1.1 


0.5851 


2.3 


3.5 








02022-2402 


TOK 41 Ba,Bb 


HIP 


9497 


9 


6709 




3 


162.2 


0.3 


0.1018 


1.0 


0.1 
















9 


6710 


R 


2 


162.7 


1.4 


0.0948 


2.9 


0.0 
















9 


7557 




3 


160.1 


4.0 


0.1007 


2.5 


0.0 








02057-2422 


WSI 71 Aa,Ab 


HIP 


9774 


9 


7557 




2 


276.1 


0.8 


0.0462 


8.2 


1.4 
















9 


7558 




2 


283.7 


2.3 


0.0553 


0.1 


1.4 
















9 


7559 




2 


287.2 


8.3 


0.0525 


1.4 


1.5 








02057-2423 


I 454 AB 


HIP 


9774 


9 


7557 




2 


154.4 


0.1 


0.8204 


1.1 


1.9 
















9 


7558 




2 


154.6 


0.3 


0.8173 


0.1 


2.2 
















9 


7559 




2 


154.7 


0.6 


0.8185 


3.1 


2.0 









Table 1 — Continued 



WDS 


Discoverer 


Other 


Epoch 


Filt 


N 


e 


P<JQ 


p 


crp 


Am 


[O-C]^ 


[0-C]p 


Reference 


(2000) 


Designation 


name 


+2000 








(mas) 


i") 


(mas) 


(mag) 




V / 


code 


02225-2349 


TOK 40 


rliJr 


1 1 ClTO 
iiU ( Z 


9.6709 


I 


2 


125.3 


1.2 


0.4969 


1.2 


3.8 
















9.6709 


R 


2 


124.8 


4.6 


0.4997 


4.6 


4.5 
















9.7558 


I 


5 


126.1 


0.9 


0.4966 


1.2 


3.6 * 








02396-1152 


FIN 312 






9.6710 


y 


2 


131.1 


0.3 


0.0762 


0.2 


0.8 


7.6 


-0.003 


Sodl999 


02442-2530 


FIN 379 Aa,Ab 


HIP 


±Z I oU 


9.6709 


I 


2 


302.2 


0.4 


0.0648 


0.7 


0.2 


-41.6 


-0.279 


Hng2005 










9.7557 


I 


3 


310.1 


0.9 


0.0727 


2.8 


0.5 


-33.9 


-0.270 


Hng2005 










9.7558 


V 


1 


312.2 


0.5 


0.0596 


0.5 


0.7 : 


-31.8 


-0.283 


Hng2005 


02460-0457 


BU 83 


rlir^ 


LZi) LZ 


9.6710 


I 


2 


14.8 


0.4 


0.9560 


0.4 


1.6 * 


1.3 


0.130 


Ole2002d 


02572-2458 


BEU 4 Ca,Cb 


"LJTT3 

rli_r 


1 Q'T^^Q 

lo / oy 


9.6709 


I 


2 


A352.4 


0.3 


0.0818 


0.9 


0.8 


10.6 


-0.985 


Sca2002c 


17031-5314 


HDS 2412 Aa,Ab 


HIP 


00401 


9.7551 


I 


4 


199.4 


0.7 


0.5939 


0.9 


2.7 * 
















9.7551 


Ha 


2 


199.3 


0.7 


0.5925 


1.3 


2.8 








17195-5004 


FIN 356 


rliJr 


HA I^Q 


9.7551 


I 


2 


A84.7 


0.2 


0.0698 


0.4 


0.0 
















9.7552 


Ha 


2 


A91.0 


0.2 


0.0647 


0.7 


0.4 : 
















9.7564 


Ha 


2 


A89.9 


2.2 


0.0645 


0.5 


0.7 : 








17248-5913 


I 385 AB 


rlir^ 


OOZ io 


9.7551 


I 


3 


121.3 


2.3 


0.3984 


4.8 


0.8 








17248-5913 


WSI 87 AD 




oOZlD 


9.7551 


I 


3 


269.4 


4.3 


0.2687 


7.7 


0.7 








18031-0811 


STF 2262 AB 


rliJr 




9.6702 


I 


2 


285.2 


0.8 


1.6053 


0.5 


0.7 * 


0.5 


-0.020 


Sodl999 










9.6702 


R 


6 


285.2 


1.2 


1.6024 


1.3 


0.7 * 


0.5 


-0.023 


Sodl999 


18112-1951 


BU 132 


HIP 


SQ1 1 A 
oy 114 


9.6702 


I 


1 


8.4 


1.3 


1.4047 


1.3 


1.3 : 








18594-1250 


KUI 89 


rli_r 


yozzo 


9.6703 


I 


2 


142.4 


0.4 


0.2159 


0.3 


0.9 


12.4 


0.019 


Msnl999c 










9.6703 


R 


2 


142.3 


0.4 


0.2173 


0.4 


1.0 


12.3 


0.020 


Msnl999c 










9.6703 


I 


8 


142.4 


0.8 


0.2163 


0.8 


0.9 


12.4 


0.019 


Msnl999c 


20311-1503 


FIN 336 


HIP 


101221 


9.6705 


I 


4 


304.7 


0.4 


0.1528 


0.6 


1.5 


-1.9 


0.000 


Msn2010 










9.6705 


R 


2 


304.8 


0.8 


0.1538 


0.6 


1.3 


-1.8 


0.001 


Msn2010 


20393-1457 


HU 200 AB 


HIP 


lu ly zo 


9.6706 


R 


2 


118.1 


0.3 


0.3313 


0.3 


0.6 


-2.1 


0.006 


Heil998 










9.6706 


I 


2 


118.1 


0.2 


0.3318 


0.2 


0.5 


-2.1 


0.006 


Heil998 


20514-0538 


STF 2729 AB 


TIJT) 

rli_r 
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* The complete list of references may be found at ^http: / /ad.usno.navy.miI/WebtextfiIes/wdsnewref.txt| 



